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We present an analysis of a scalar field model of dark energy with an exponential potential using 
the Dark Energy Task Force (DETF) simulated data models. Using Markov Chain Monte Carlo 
sampling techniques we examine the ability of each simulated data set to constrain the parameter 
space of the exponential potential for data sets based on a cosmological constant and a specific 
exponential scalar field model. We compare our results with the constraining power calculated 
by the DETF using their u wq — w a " parametrization of the dark energy. We find that respective 
increases in constraining power from one stage to the next produced by our analysis give results 
consistent with DETF results. To further investigate the potential impact of future experiments, 
we also generate simulated data for an exponential model background cosmology which can not be 
distinguished from a cosmological constant at DETF "Stage 2", and show that for this cosmology 
good DETF Stage 4 data would exclude a cosmological constant by better than 3a. 

PACS numbers: 



I. INTRODUCTION 

In the late 90 's, two independent teams presented evi- 
dence from supernova observations that the universe, in- 
stead of slowing down due to gravity, is accelerating [TJH]. 
In the standard cosmological framework, the acceleration 
is caused by a mysterious new form of matter, dubbed 
"Dark Energy" , that makes up roughly 70% of the uni- 
verse. There is a wide variety of possible explanations for 
dark energy. The simplest model that provides a good 
fit to the data is a cosmological constant. A cosmologi- 
cal constant is equivalent to a homogeneous fluid with a 
constant energy density and a ratio of pressure to energy 
density (the "equation of state parameter" w), equal to 
— 1 at all times. Yet, despite compelling evidence for the 
existence of dark energy, it is unclear whether the dark 
energy density is constant or varies with time. There are 
many different proposals for a dynamical form of dark 
energy, one of them being quintessence. Quintessence 
describes the acceleration being caused by a scalar field, 
</>, but even just among quintessence models there is a 
tremendous variety of possible behaviors. There is con- 
siderable interest in acquiring better data in order to im- 
prove our understanding of dark energy. 

Recently the Dark Energy Task Force (DETF) re- 
leased a report charting a course for future experiments 
[3J. They modeled dark energy as a homogeneous and 
isotropic fluid with an equation of state parameterized by 
w(a) = wq + w a (l — a), where the scale factor a = 1 to- 
day. Defining Stage 1 to be what is already known, they 
forecasted data for three additional experimental stages: 
Stage 2 data represents on-going experiments that will 
be completed in the near future. Stage 3 data sets repre- 
sent medium sized proposed experiments. Lastly, Stage 4 
data sets represent proposed large scale future space and 
ground-based experiments. Each stage is further catego- 
rized as either "optimistic" or "pessimistic" depending on 
how well the systematics are expected to be constrained. 



The scientific impact of a stage was quantified in terms 
of a "Figure of Merit" (FoM). The Figure of Merit is de- 
fined to be the ratio of the area of the 2a contour of the 
wo — Wa space for Stage 2 divided by the area of the 2a 
contour of the wq — w a space for Stage 3 (or 4) . 

The DETF analysis leaves several open questions, 
some of which our research seeks to address. The wq — w a 
parameterization is not motivated by a physical model of 
dark energy and provides cosmological solutions that may 
be very different from a scalar field model. As illustrated 
in Fig. [TJ the w curves generated by the exponential 
scalar field model that we consider in this paper are not 
especially well fit by curves in the wq — w a family, except 
for those nearly identical to w = — 1. Thus the relation- 
ship between the DETF results and the impact of future 
experiments on scalar field models is not clear. A good 
way to clarify this point is to model the impact of future 
data sets directly on particular scalar field quintessence 
models, which is what we do here. Our work comple- 
ments the DETF report as well as work by other authors 
using alternative w(a) parameterizations [HIS], parame- 
terizations oi pde [3], and model independent scalar field 
parameters [7]. 

The exponential scalar field model has been used in 
many different cosmological contexts due to its ability to 
give scaling solutions for the scalar field energy density 
P0 where d ^og{a)) ~^ ^ e cons t an t depends on pa- 
rameters in the scalar field potential as well as the other 
forms of matter present in the universe . Originally the 
potential was used for power law inflation models and was 
shown to have a range of attractor solutions |8J . Its abil- 
ity to produce attractor solutions that scale like the back- 
ground energy density made it an interesting choice for a 
dark matter candidate [SJ [ID1 |TT] . The variety of scaling 
solutions is well covered by Copeland, et al. |T2j where 
they argue that consideration of "fine tuning" parameters 
and constraints from nucleosynthesis give A > 20 as a 
natural choice. This range of A values wouldn't allow for 
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FIG. 1: Here we illustrate the differences between w(a) curves from the exponential model (solid) and from the ansatz used 
by the DETF (dot-dashed). Curves are plotted with respect to a (left panel) and z (right panel). The four curves from top 
to bottom are given by (w ,w a ): (-0.77673,-0.2327), (-0.9211,-0.0789), (-0.9797,-0.0203), (-0.9992,-0.0008). Using 
the final parameterization (explained in Sect. Ill ) , the exponential model parameters are given by (A, Vi): (0.07,0.3725), 
(0.35,0.42), (0.7,0.42), (1.2,0.52) (in units defined in the text). Both models give w(a — 1) to be the same value for each set 
of curves and w{ai) = — 1 for all curves shown. 



late time cosmological acceleration and is therefore ruled 
out as an explanation for dark energy. The "fine tuning" 
that is required to successfully describe dark energy with 
this model is needed so that the scalar field energy den- 
sity can be initially very small; of the order of the dark 
energy density today. This has caused this model to be 
discarded by many authors on the basis that the model 
has lost the theoretical generality that made the poten- 
tial initially interesting. However, as a practical matter 
the fine tuning is straightforward to implement, and the 
potential is very simple and easy to work with. Because 
of this simplicity, we found it valuable to have this poten- 
tial as part of our larger project (which includes a vari- 
ety of more complicated quintessence potentials [13j[T4]). 
The simplicity helped us deal with a number of technical 
issues first with the exponential model and then trans- 
fer our understanding to the more complicated cases. In 
addition, realistic cosmologies for the exponential model 
have their special forms for w(a) (illustrated in Fig. [T]). 
We found it useful to include this family of w(a) curves 
in our set of possibilities to more fully understand the 
constraining power of future data sets. 

The paper is organized as follows. In Section [IT] we 
provide an introduction to our scalar field model and its 
cosmological solutions. In Section [ITT| we describe how we 
come about our choice of parameterization, as this is a 
critical step in the MCMC analysis. (An account of our 
general MCMC methods and data modeling can be found 
in the appendix of our companion paper |13J . This pa- 
per contains only information specific to the exponential 
model.) Section IV presents our results for data simu- 



lated using a background cosmology with a cosmological 
constant and then Section [V] presents results where the 
data is based on a cosmology with exponential model 
quintessence. Finally, we summarize our key results in 
the conclusions. 



II. EXPONENTIAL MODEL COSMOLOGY 

We model dark energy as a homogeneous scalar field 
evolving in an exponential potential 



(1) 



The cosmological evolution of this scalar field in a FRW 
universe is then given by solving: 
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where M p is the reduced Planck mass. The equation of 
state of the scalar field is given by 
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FIG. 2: The left panel shows three examples of the exponential potential (dashed line 
long-short dashed line: A_^0.35,!/ = 0.32,0/ = -0.78, solid line: A = 0.7, V = 0.38,0/ 

The path of the scalar field is depicted by thick solid curves 



A = 0.07, Vo = 0.24,0/ = -6.28, 
= —0.14. Vq is given in units of h 2 , 
The corresponding w behavior is 



as mentioned in Section IV I 
shown in the right panel. The solid curve gives the potential used for the fiducial model discussed in Section [V] 



which we will use in discussing an evolving dark energy. 
In this picture a cosmological constant is equivalent to a 
scalar field with w = — 1. 

In our analysis we initially set ^ = |18j . This leaves 
the dynamics of the field completely determined by the 
slope and curvature of the potential: 

(6) 

Since the initial field velocity is zero, the initial equation 
of state is w = —1, mimicking a cosmological constant. 
As the field begins to roll the equation of state begins 
to depart from —1. The rate of this departure is deter- 
mined by the steepness of the potential. A steeper slope 
gives changes in <fi that correspond to larger changes in w. 
Likewise, a flat slope gives little change in <fi and there- 
fore has a cosmology similar to a cosmological constant, 
as shown in Figure [2| 

The possible scaling solutions achievable by the expo- 
nential model are systematically discussed by Copeland, 
et al. 12J. If a scaling solution is reached before the 
onset of dark energy domination the universe will not 
accelerate and therefore is a poor match to current data. 
However, there are subsets of scaling solutions that reach 
their scaling solution after dark energy domination and 
provide different fates for the universe. These scaling so- 
lutions fall into two categories: (i) < A < \/2 or (ii) 
V2 < A < A*. Solutions with A values in category (i) 

approach scaling where w — > 4j — 1, giving late time ac- 
celeration. For category (ii), we define A* to be the value 
that gives w(a = 1) < — |, but the scaling solution leads 
to w(a > 1) > — |. It is possible for A* to be larger than 
V3 and therefore have w(a > 1) — > 0. This value depends 



on the initial scalar field energy density, p^j, which in 
turn determines when the field begins to approach its 
scaling solution, i.e. when the field starts rolling. We are 
allowing p^j and other cosmological parameters to float 
so a universal value of A* cannot be uniquely determined. 



III. PARAMETRIZATION 

In order to run our MCMC analysis on the potential 
V = Vbe _A ^, we first need to make a careful choice of 
parameters. The obvious choice of the potential param- 
eters Vo, A, and the initial field value </>/ presents several 
problems. Rewriting the potential V = Vbe" A ^ — ► V = 
e ln(v o )-A0 reveals a degeneracy between ln(Vb) and \4>i- 
For fixed values of A, a change in <pj and a corresponding 
change in ln(Vb) gives identical cosmological solutions. 
This will lead to an unconstrained and uninteresting pa- 
rameter space. Fixing Vq removes this degeneracy. We 
make the choice of V — pa = 8.74 x 10~ 121 (in reduced 
Planck units), which is the value of the cosmological con- 
stant energy density used by the DETF. This is the sim- 
plest choice, although not absolutely necessary. Other 
choices of Vq would provide equivalent cosmological so- 
lutions. 

Removing the degeneracy and fixing Vo leaves A and 
4>i as the two model parameters. However, this choice 
leads to an "infinite direction" in A — 4>i space: Since the 
data for the first part of our analysis is modeled on a 
cosmological constant, the most probable values of A are 
those where A approaches zero. As A approaches zero, 
cj>i can take any value and produce solutions indistin- 
guishable from a cosmological constant. This leads to an 
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FIG. 3: Likelihood contours in the Vi — A space for cosmological constant data models. The three contours give 68.27%, 95.44%, 
and 99.73% confidence regions. 



infinite unconstrained direction in parameter space that 
is uninteresting and also fatal to the MCMC techniques. 

One can resolve this problem by placing a bound on A 
or (f)j. For small values of A, a bound placed on A is nearly 
equivalent to placing a bound on w(a =1). A choice 
of a bound on A can be chosen such that the difference 
from w = — 1 is small, however the choice is arbitrary. 
Further, for data based on a A universe, the closer the 
bound is placed to A = 0, the more the allowed region of 
parameter space squeezes against this bound as smaller 
values of A allow a wider range of (f>j, basically partially 
restoring the degeneracy we are trying to eliminate. This 
arbitrariness and distortion of allowed parameter ranges 
make bounding A a poor choice for addressing the pa- 
rameter space degeneracies in this model. The squeezing 
effect leads to incorrect conclusions about allowed values 
of A. The space appears to disfavor larger values of A, 
or equivalently larger departures from w(a = I) = —I, 
than would the space in the final parameterization that 
we discuss next. 

We find the best choice of free model parameters to be 



Vj, where Vj = V(4>i), and A. The value of <pj is then 
determined from A and Vj. This parameterization avoids 
the degeneracy discussed in the previous paragraph since 
a cosmological constant of a particular value is only rep- 
resented at one point (Vj = Pa and A = 0) in the A — V/ 
space. Values similar to a cosmological constant are ex- 
plored without an arbitrary bound placed on any param- 
eter. This allows the MCMC method freedom to explore 
a more natural space. 

There is no loss of generality with this choice of param- 
eterization as is easily seen by the slope and curvature of 
the potential in the new parameters: 



dV 
~d4 



-AV/ 



d 2 V 
d<t> 2 



A 2 Vj 



(7) 



This parameterization allows a simple intuition about the 
role of these parameters forming the solutions. Small 
values of A give a flat potential and the scalar field will 
be stationary, independent of the choice of V/. For large 
values of A the field will roll and the amount to which it 
does will depend on the value of V[. 
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FIG. 4: Likelihood contours in the A — 8lude- Deviations of Suide from indicate evolving dark energy. The background 
cosmology has a cosmological constant. The boxes in the top left panel (Stage 2) respectively show the size of the axes for 
Stage 3 and Stage 4 plots in Fig. [5] The contours give 68.27%, 95.44%, and 99.73% confidence regions. 



IV. COSMOLOGICAL CONSTANT FIDUCIAL 
DATA 

With our parameterization firmly in hand, we now an- 
alyze the DETF data sets based on a cosmological con- 
stant cosmology using the MCMC technique. The like- 
lihood contours for Stage 2, Stage 3 Photo Optimistic, 
Stage 4 Ground LST Optimistic, and Stage 4 Space Op- 
timistic are shown in Fig. [3] for the A — Vj space. In all 
plots shown in this paper we use the DETF supernova, 
weak lensing, baryon oscillation, and PLANCK (using 
the alternate parameters as in 0) data sets but not the 
cluster data sets due to technical problems adapting the 
DETF cluster data models to our methods. These tech- 
nical problems are similar to those outlined in [5 , where 
similar issues where encountered. Our plots were con- 
structed by marginalizing over all the cosmological pa- 
rameters, ui m , tuk, lob, #0 n S' t ne various nuisance pa- 
rameters, and/or the photo-z parameters. The nuisance 
and photo-z parameters are detailed in the appendix of 



one of our companion papers [13J. The fiducial values for 
the cosmological parameters are shown in Table [TJ The 
values for all energy densities and Vj in the remainder of 
the paper are in units of h 2 , where h = 

Figure [4] gives likelihood contours in A — 6oj ee space, 
where 5lo de = u DE (a = l)-uj DE {ai). Here lo de = j^h 2 
and p c — 3M 2 H 2 . The value of Sujde gives the amount 
the dark energy density has changed since the simulation 
started at scale factor aj (well in the radiation era) . Val- 
ues of Slu y) e different from zero correspond to dynamical 
dark energy. Figure [5] shows an enlarged version of the 
A — 5u>d E space for Stage 3 and Stage 4 experiments. 

For Stage 2, values of A from to about 0.15 in Fig. 
[4] correspond to shallow slopes and don't allow for much 
change in the amount of dark energy. For these values 
of A there is a spread in Vj in the A — Vj space (Fig. |3j. 
Since these values are consistent with a non-evolving dark 
energy, the spread in Vj is essentially a measure of how 
well the experiment is measuring L0 EE (a =1). Values of 
A > 0.15 all correspond to detectable differences from a 
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FIG. 5: Enlarged plots of the A — Suode confidence contours for cosmological constant data models. The three contours 
correspond to 68.27%, 95.44%, and 99.73% confidence regions. The axes correspond to the boxes in Fig. [4] 



cosmological constant. This portion of the A — Vj space 
has an upturned feature. As the slope gets steeper the 
field needs to start higher up in the potential in order 
to roll down to acceptable values of u>de{<i = 1)- These 
features are evident in the plots for Stage 3 and Stage 4, 
although by Stage 4 it is less clear as the parameter space 
has shrunken to values closer to a cosmological constant, 
and so the upturned trend has diminished. 

Comparing Stage 2 to Stage 3 Photo Optimistic, and 
then on to Stage 4, there is a significant tightening of 
the allowed area in parameter space. This increased con- 
straining power is similar to the factors of about three 
(Stage 2 to Stage 3) and ten (Stage 2 to Stage 4) in- 
crease in constraining power noted by the DETF in the 
wo — w a space. The A — Suj^e contours in Fig. [4] al- 
low one to interpret the constraining power already seen 
in the A — Vj space in terms of a specific aspect of the 
dark energy dynamics, namely the overall change in dark 
energy density (given by 6u>de)- 

Figure [6] shows plots of the allowed functions of w(z) 



for each data set. The plots are constructed by selecting 
around 100 points taken uniformly within the 3er contour 
(thus also including points within the 1 and 2a contours). 
This is done to illustrate the full range of solutions not 
excluded at better than 3er. The furthest most curve 
from w = — 1 for each Stage corresponds to the top right 
most tip in the A — Vj space. 

Since we are interpreting the data using the exponen- 
tial quintessence model, we can use our knowledge of how 
the cosmological solutions vary with A to discuss con- 
straints on future cosmology as well. Since all the stages 
depicted in Figure [3] strongly favor A < \/2, essentially 
all future w(a) behavior consistent with these data sets 

A 2 1 

will approach scaling solutions where w — > -3 — 1 < — | 
and give accelerating universes where — » 1. The scal- 
ing solutions have not been reached by today but will be 
approached in the future. Therefore, by Stage 2, solu- 
tions that lead to universes with a non-accelerating fate 
have been ruled out in this scenario. We will revisit this 
point in the next section in the context of a different 
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FIG. 6: The w(z) behavior for a sample of points covering the 
constant. 

background cosmology. 



V. EXPONENTIAL MODEL FIDUCIAL DATA 

Next we consider the case where the universe happens 
to have dark energy described by the exponential model. 
We select a particular fiducial model of dark energy to 
illustrate the potential impact of Stage 4 data. We use 
exponential model parameters of A = 0.7 and Vj = 0.42 
for the fiducial model (the remaining parameters were 
given the same values we use for the cosmological con- 
stant background model, given in Table 0. 

We choose the fiducial values by finding a point in the 
X — Vi space in Fig. |3]for Stage 2 that was just outside a 
lcr detection but was excluded by better than 3er in Stage 
4 Optimistic (for both ground and space). This point 
corresponds to w(a = 1) = —0.92, as shown in Fig. [2] 
The results of Stage 2, Stage 3 Photo Optimistic, Stage 4 
Space Optimistic and Stage 4 LST Optimistic are shown 
in Figure [7] for the A — Vi space. The A — Slode space is 
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shown in Fig. [8] Fig. [9] shows an enlarged version of the 
A — Slude space for Stage 4 experiments. 

For Stage 2, the likelihood contours in X—Vi space look 
similar to the case with a cosmological constant back- 
ground cosmology. The range of A is nearly the same 
with all A < Therefore the same conclusions about 
strongly favoring future scaling behavior that were made 
in the previous section apply here. The upturned trend 
is a little more dramatic since dark energy solutions with 
more evolution are favored. The major difference in the 
space is that values of A < 0.05 are outside the la con- 
tour. Comparing this with the la region of 5lude shows 
that these values are consistent with a non-dynamical 
dark energy. It is not until we reach values of A f=a 0.2 
that we find a la region that corresponds to an evolving 
dark energy in the Sujde plot. Therefore, even though 
the parameters consistent with a cosmological constant 
fall outside the la contour in the A — Vj plot, a non- 
dynamical dark energy, and thus a cosmological constant, 
is not ruled out at even the la level by analyzing the con- 
tours in the A — Swde space. This apparent discrepancy 
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FIG. 7: Plots of the likelihood in Vi — A space for data based on the fiducial exponential background cosmology. The three 
contours correspond to 68.27%, 95.44%, and 99.73% confidence regions. 



between the two pictures is in fact a common situation 
when examining relatively low likelihood contours in dif- 
ferent parameter spaces. More clear signals will only be 
found when looking at phenomena that are rejected at 
higher likelihood levels. 



TABLE I: ACDM (left column) and Exponential (right col- 
umn) Fiducial Parameter Values. The nuisance and photo-z 
parameters are 0. 



ujde 0.3796 


0.3796 


co m 0.146 


0.146 


LO k 0.0 


0.0 


ujb 0.024 


0.024 


n s 1.0 


1.0 


8 C 0.87 


0.87 


A 0.0 


0.7 


Vi 0.38 


0.42 



By Stage 3 Photo Optimistic, the cosmological con- 
stant is now excluded in the la contour although well 
within the 2a contour in the A — Vi space. This is con- 



sistent with the A — Sojde plot in Fig. [8| The increased 
constraining power is again equivalent to the DETF re- 
sult for the X — Vi space. However, the range of A has not 
changed much within the 3a contour, allowing the range 
of evolving dark energy solutions to be nearly as large as 
Stage 2, as seen in Fig. [8] 

Stage 4 clearly differentiates between the exponential 
fiducial model and a cosmological constant by better than 
3ct, as shown in both X — Slvde and A— Vi spaces depicted 
in Figures [7] and [8] Again, the increased constraining 
power is consistent with the DETF results for both Stage 
4 experiments. 



Figure 10 shows the span of the w(z) solutions for each 
of the Stages determined in the same way as in the pre- 
vious section. The curve with the greatest w(z = 0) 
departure from —1 is obtained from the top right tip of 
the 3a curve in the X — Vi space. For the Stage 4 plots, 
the bottom most curve, corresponding to the closest ap- 
proach to a cosmological constant, corresponds to the 
bottom left tip of the 3a region in the A — Vi space. 
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FIG. 8: Plots of the A — 5u>de confidence contours for data based on the fiducial exponential background cosmology. The three 
contours correspond to 68.27%, 95.44%, and 99.73% confidence regions. The box in the top left figure (Stage 2) shows the size 
of the axes for the Stage 4 plots in Fig. [9] 



VI. DISCUSSION AND CONCLUSIONS 

We have analyzed the exponential scalar field model 
using MCMC techniques for the DETF simulated data 
sets representing future dark energy experiments. We 
have demonstrated the ability of these experiments to 
place significant constraints on the parameters of a scalar 
field model. The relative constraints on the size of the 
A — Vi space between various data sets produce values 
similar to the constraints computed by the DETF in the 
wo — w a space. In addition to placing constraints on 
the quintessence parameters, we also presented our re- 
sults in terms of the evolution of the dark energy in our 
quintessence model. This allows us to distinguish more 
directly between a cosmological constant and our partic- 
ular scalar field model. We have presented plots for a 
characteristic selection of combined DETF data models, 
but in the course of this work we have also examined 
similar plots for a much wider range of DETF data mod- 
els, including data models representing single techniques. 



We found that the consistency with constraints reported 
by the DETF in wq — w a space to hold across the entire 
range of data choices and combinations we considered. 

We based our data on two different background cos- 
mologies, one with a cosmological constant and one with 
exponential model quintessence with specific parame- 
ters. We found that the equivalence with the DETF 
results held in both cases. Our specific background 
quintessence model was chosen (with the parameter val- 
ues of V(4>i) = 0.42 and A = 0.7) in order to illustrate 
the power of Stage 4 experiments. For this model, the 
maximum deviation from w = — 1 occurs today, with 
w(a — 1) = —.92. We found that if the universe is ac- 
celerating due to this particular exponential quintessence 
model then a cosmological constant dark energy model 
can be ruled out to at least 3a by good Stage 4 experi- 
ments. For this background cosmology, the cosmological 
constant is within the la contour at Stage 2 and the 2a 
contour at Stage 3. 

We note that there are a number of ways experiments 
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FIG. 9: Enlarged plots of the A — 5lude confidence contours for data based on an exponential model. The three contours 
correspond to 68.27%, 95.44%, and 99.73% confidence regions. The axes correspond to the boxes in Fig. [8] 



might be optimized to do better than the cases consid- 
ered by the DETF (see for example [THE)]). We have not 
included such ideas in our work, with an eye for offering 
more direct points of comparison with the DETF. How- 
ever, improvements such as these could lead to more pow- 
erful constraints on quintessence models than we have 
calculated here. 

We have found in this work and in our companion pa- 
pers |13l I14) that a wide variety of quintessence mod- 
els (with widely varying families of functions w{z)) are 
constrained by DETF data in a way comparable to the 
constraints found in wo — w a space by the DETF. As 
discussed in [17], we believe that this is related to re- 
cent work by one of us (AA) and Bernstein [5], where 
it was demonstrated that, overall, the good DETF sim- 
ulated data sets could constrain significantly more than 
two dark energy parameters. From this point of view 
our various models of dark energy are just sampling dif- 
ferent more or less "random" combinations of the "well 



measured modes" discussed in [5] and in each case are 
coming up with similar results. 
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